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PERFORMANCE  OF  THE  OPTIMUM  AND  SEVERAL  SUBOPTIMUM 
RECEIVERS  FOR  THRESHOLD  DETECTION  OF  KNOWN  SIGNALS 
IN  ADDITIVE,  WHITE,  NON-GADSSIAN  NOISE 


INTRODUCTION 


The  additive  noise  encountered  at  a  receiver  input  is  often  non-Gaussian. 
If  the  non-Gaussian  nature  of  the  noise  is  lot  taken  into  account  in  the 
design  of  the  receiver,  significant  performance  degradation  of  the  receiver 
can  be  expected.  The  purpose  of  this  technical  report  is  to  describe  the  per¬ 
formance  that  can  be  expected  from  the  optimum  and  several  subopt inua  receivers 
used  for  detecting  known  threshold  signals  in  an  additive,  white  (i.e., 
statistically  independent  noise  samples),  non-Gaussian  noise  environment. 

The  optimum  receiver  structure  for  detecting  known  threshold  signals 
in  additive,  white,  non-Gaussian  noise  is  the  same  as  that  which  should  be 
used  if  the  noise  were  Gaussian,  except  that  a  zero  memory  nonlinearity  is 
placed  between  the  receiver  input  and  the  Gaussian  detector.  The  input- 
output  characteristic  of  the  nonlinearity  is  given  by  -d/dx  (In  pn(x)L  where 
Pn(x)  is  the  density  function  of  the  noise  alone.  The  measure  of  performance 
improvement  obtained  by  including  the  nonlinear  device  is  given  by  the  ratio 
of  the  signal  to  noise  ratio  (SNR)  at  the  receiver  output  with  the  nonlinear¬ 
ity  in  the  circuit  to  the  SNR  at  the  receiver  output  without  the  nonline¬ 
arity.  The  magnitude  of  this  improvement  is  evaluated  here  using  Middleton's 
Class  B  Noise  Model.  The  performance  of  several  more  easily  implemented 
subopt imum  nonlinearities,  i.e.,  clipper,  hole  runcher,  and  hard  limiter, 
are  also  evaluated  using  Middleton's  model;  and  typical  performances 
relative  to  the  optimum  are  presented.  Finally,  the  amplitude  probability 
distributions  (APD's)  of  high  level  extremely  low  frequency  (ELF)  noise  are 
given  and  compared  with  the  APD's  resulting  from  Middleton's  model.  From 
these  APD's  and  the  parametric  results  presented,  the  expected  performance 
of  the  optimum  and  subopt iiaum  nonlinearities  is  derived. 


OPTIMUM  RECEIVER  STRUCTURE 


The  optimum  receiver  structure  for  the  threshold  detection  of  binary 
coherent  signals  in  an  additive,  white-noise  environment  is  well  known.1  3 
The  derivation  of  this  structure  is  summarized  here. 

Consider  a  discrete  time,  memoryless  communication  channel  in  which  one 
of  two  known  signals,  Si  =  ***,  S^] '  or  S2  =  [S21,  S22*  ***,  S2n)  r, 

is  transmitted  every  T  seconds  by  successively  transmitting  the  signal  compo¬ 
nents,  Sji,  (i  =  1,  2,  — ,  N),  every  At  seconds.  Each  digit,  zj,  in  the 
received  sequence  of  digits  z,  is  the  sum  of  one  of  two  known  signals  and  a 
white-noise  component,  ni.  The  receiver,  then,  must  decide  between  the  two 
hypotheses.  Hi:  _z  =  Si  +  n^  and  H2:  £  =  §2  +  n,  in  such  a  manner  that  some 
criterion  is  optimized.  For  example,  one  may  wish  to  minimize  the  expected 
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risk  (Bayes  Test)  if  costs  have  been  assigned  to  all  joint  outcomes  of 
correctly  or  incorrectly  choosing  a  particular  hypothesis.  Alternatively, 
lacking  such  a  cost  structure,  one  nay  minimize  the  probability  of  an  in¬ 
correct  choice  (Ideal  Observer  Scheme). 

Because  the  noise  is  white  and  the  signal  components  are  known  com¬ 
pletely,  the  conditional  probability  density  function  of  the  received  vector, 
z ,  can  be  expressed  in  terms  of  the  density  functions  of  the  individual  noise 
components,  pR(a^);  i.e., 

N 

P«<*/V  =  ?z(zr  V  *"*  W  =  “  Pn(zi  -  sji3  (j  =  2>2)  * 

—  —  j  i=i 


Hie  receiver  decides  between  and  H2  by  comparing  pz (z/Hj)  with 
pz(z/H2)  relative  to  some  constant  threshold.  A,  determined  by  the  a  priori 
probabilities  of  Hj  and  H2  and  the  optimality  criterion.  This  is  to  say  that 
the  receiver  forms  the  following  likelihood  ratio  as  its  test  statistic: 


n  pn(^  -  su)  h 

i=l  n  1  11  >A  x 

jj,  »„lzi  -  s2i>  "2 


Since  the  logarithm  function  is  monotonic,  an  equivalent  test  is 


Considering  the  case  of  threshold  detection,  let  the  signal  components, 
Sjj(j  =  1,2,  i  =  1,2  — ,  N),  be  sufficiently  ssall  so  that  each  term  in  the 
above  expression  can  be  expanded  in  a  first  order  Taylor  series  about  z.  . 
The  log  likelihood  ratio  test  then  becomes  1 


In  A  = 


N 

T|  Ir.  p  (z.  -  S_.)  -  V  In  p  (z.  -  S„.) 
f-,  1  lx  J-i  rn  1 
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i=l 


i=l 


$ 
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In  A 


In  A 


N  r 
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S2i 


dz. 
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ln  A  , 
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or 


In  A  =  £ 


i=l 


»  N  m  V 

ln  Pn'V)  -  £  s2i(-  ar ln  PnC2i’)  j1 

*  H- 


In  X  .  (1) 


Equation  (1)  implies  the  block  diagram  of  figure  1  for  the  optimal 
receiver  structure.  Note  that  it  is  just  a  correlation  receiver  in  which 
a  zero  memory  function,  -  d/dzj  (ln  Pn(zi))>  operates  on  the  received  signal 
prior  to  the  correlation.  Note  also  that  the  receiver  is  canonical;  i.e., 
its  general  fora  and  decision  logic  is  not  based  on  any  particular  model  of 
the  noise  environment.  However,  the  receiver  must  "know"  (or  estimate)  the 
noise  density  function  in  order  to  implement  the  -  d/dzj  (ln  Pn(zi)  operation 
in  the  block  diagram.  If  the  noise  environment  is  changing  slowly  in  time, 
the  optimal  receiver  must  be  adaptive. 


When  the  noise  is  Gaussian,  the-d/dZj(ln  pn(zj))  function  is  linear  in 
z j  and,  as  such,  may  be  incorporated  as  a  gain  factor  (weight)  within  the 
signal  component,  Sjj,  so  that  the  optimum  receiver  structure  is  exactly 
that  of  a  correlator: 


Gaussian  noise 


m  = 


Z. 

1 


For  a  non-Gaussian  noise  density,  however,  the  derivative  of  the  log  density 
function  is  not  a  linear  function  of  the  received  waveforms.  Thus,  the 
optimum  receiver  structure  is  a  correlator  preceded  by  a  zero  memory  circuit 
whose  input-output  characteristics  depend  only  on  the  noise  probability 
density  function. 


As  an  example,  it  has  been  proposed  that  naturally  occurring  electro¬ 
magnetic  noise  below  100  MHz  is  Gaussian  distributed  with  occasional  large 
amplitude  noise  bursts,  which  occur  often  enough  to  affect  the  tails  of  the 


noise  densil Several  authors  6  have  suggested  that  the  density  function 
describing  t  ese  tails  follows  an  e-aiz!  (0  <  k  <  1)  law.  The  smaller  the 
k,  the  more  impulsive  the  noise  (see  figure  2a).  For  this  noise  density, 
the  input -output  characteristic  of  the  optimal  nonlinearity  is  shown  in 
figure  2b.  As  the  noise  becomes  more  impulsive,  the  nonlinear  device 
suppresses  the  larger  noise  excursions  to  a  greater  degree. 


PERFORMANCE  MEASURE 


After  the  optimal  receiver  structure  has  been  derived,  a  natural  question 
to  ask  is  how  much  performance  improvement  can  be  expected  using  this  receiver 
compared  with  one  in  which  the  nonlinearity  prior  to  the  correlator  is  sub¬ 
opt  iaal  or  absent?  Usually  the  performance  of  a  receiver  is  measured  by  the 
probability  of  error  for  a  given  signal  energy.  This  requires  that  the 
probability  density  function  of  the  test  statistic  be  known  or,  alternatively. 
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that  the  probability  of  error  can  be  bounded  in  some  manner.  Under  the 
assumption  of  threshold-type  signals,  the  correlator  must  sum  many  received 
symbols  to  increase  the  SNR  for  satisfactory  performance.  As  was  seen 
previously,  the  nonlinear  device  tends  to  suppress  large  amplitude  vari¬ 
ations;  thus,  each  received  digit  is  bounded.  (From  a  practical  point  of 
view,  the  received  symbols  are  always  bounded  owing  to  the  finite  responses 
of  amplifiers,  filters,  etc.,  even  when  the  receiver  is  linear.)  The  test 
statistic,  therefore,  is  the  sum  of  independent,  finite  variance  random 
variables  and,  by  the  Central  Limit  Theorem,  is  asymptotically  Gaussian 
distributed.  Since  the  density  function  of  a  normal  random  variable  is 
completely  specified  by  its  mean  and  variance,  the  SNR  at  the  output  of  the 
correlator  may  be  used  to  express  the  receiver  performance.  Hence,  define 
the  improvement  factor,  I,  as  the  ratio  of  the  SNR*s  at  the  output  of  the 
receivers  with  and  without  a  nonlinear  noise  processor: 


snrnl 


SNR  at  output  of  receiver  with  nonlinearity  present 


SNR  at  output  of  receiver  with  nonlinearity  removed 


Assuming  threshold  binary  equal-energy  signals,  the  improvement  factor 
for  the  optimal  receiver  is  (see  appendix  A) 


-v  4 

I  =  a 
opt 


r  r  d ,  ,  ,12  , , ,  2r[^pnf2)] 

L  [-  ai ln  f„wj  Pnw  dz  * a  L. — rnr 


dz  ,  (2a) 


where  c“  is  the  variance  of  a  received  symbol  (i.e, ,  the  power  of  the  noise 
component).  When  the  noise  is  Gaussian,  the  improvement  factor  equals  one. 

Although  the  optimal  receiver  structure  is  canonical  in  form,  it  may  be 
difficult  to  implement.  This  leads  one  to  investigate  simpler,  suboptinal 
nonlinear  processors  using  the  improvement  factor  to  compare  the  performance 
of  various  processors.  Appendix  A  derives  the  improvement  factor  for  three 
common  nonlinear  devices:  the  hard  limiter,  the  clipper,  and  the  hole 
puncher.  The  results  are  listed  here  for  convenience: 


*HL  *  4°2  P»™ 


H 


I c  =  p2o2y/ 12  z2  Pn(z)  dz  +  c2(l  -  p)j 

IhP  "  °2[P  *  2c  PnH2/[2  Jo  *2  Pn(2)  d2]  * 


4 


■wr-ssaESl 
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where 


p  =  |  p  (z)  dz  =  Probability  of  being  in  linear  range  of  nonlinearity 
-c 

c  =  Linear  range  of  nonlinearity  (see  figure  3). 

The  input-output  characteristics  of  the  above  devices  are  shown  in 
figure  3.  Note  that  the  hard  limiter  is  the  optimum  nonlinearity  when  the 
noise  has  a  density  function  of  the  form  Pn(z)  =  Ae~ i z  L  Comparing  figure  2b 
with  figure  3,  one  sees  that  the  clipper  is  optimal  and  that  the  hole  puncher 
is  near  optimal  when  the  noise  is  Gaussian  with  tails  distributed  as  e-Jzl 
for  k  =  1  and  0,  respectively.  It  should  also  be  pointed  out  that  although 
the  noise  power,  a2  ,  appears  in  each  of  the  above  expressions,  it  does  not 
affect  the  comparison  of  one  nonlinear  device  with  another,  though  it  cer¬ 
tainly  is  important  when  comparing  a  nonlinear  receiver  with  a  linear  one. 


CANONICAL  NOISE  MODEL 


As  is  evident  in  the  above  discussion,  the  optimal  receiver  structure 
is  highly  dependent  on  the  probability  density  function  of  the  noise. 

Although  various  noise  models  have  been  proposed  for  particular  environments, 
the  only  general  model  available  to  date  is  that  suggested  by  Middleton.7  9 
His  model,  the  Class  B  Noise  Model,  which  is  analytically  tractable  and  is  in 
excellent  agreement  with  measured  data  for  a  variety  of  cases,  is  canonical 
in  nature.  That  is  to  say  the  noise  density  function,  which  is  characterized 
by  various  parameters,  does  not  change  form. 

Based  on  the  bandwidth  of  the  interference  relative  to  the  receiver 
bandwidth,  Middleton  defines  two  major  classes  of  noise.  Of  interest  here 
is  the  Class  B  Noise  Model  in  which  the  interference  is  highly  impulsive, 
resulting  in  a  noise  bandwidth  much  greater  than  the  receiver’s  front  end 
bandwidth.  For  this  case,  the  model  for  the  amplitude  noise  probability 
density  is  expressable  as  two  distinct  functions:  one  valid  for  |z|  <  zQg 
and  the  other  for  |zj  >  z0g  .  The  second  function  approaches  zero  much 
more  rapidly  than  the  first  as  z  +  ®  .  It  was  already  pointed  out  that  a 
nonlinear  device  operating  on  threshold  signals  in  this  type  of  noise  envi¬ 
ronment  will  suppress  large  voltage  excursions  in  the  received  waveform. 

Thus,  the  noise  power  at  the  output  of  the  nonlinear  device  is  primarily 
determined  by  the  noise  density  function,  Pn(z),  for  |z|  less  than  zc  , 
i.e.,  the  value  at  Which  significant  suppression  occurs.  The  more  impul¬ 
sive  the  noise  environment,  the  more  likely  zc  is  much  smaller  than  z0j$  . 

For  this  reason,  only  the  first  function  in  Middleton’s  model  is  used  here. 
Thus, 


P„W 


=  PrU) 


si?0 


ma  +  1 


1 

2 


-*) 


|z|  ^ 


'oB 


(3) 
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where 

0  <  a  <  2.,  A  >0 
a 

T(*)  =  Gamma  Function 

jFj(a;b;x)  =  Confluent  Hypergeometric  Function. 

Three  values,  Aa  ,  a  ,  and  zqB  ,  define  the  noise  model.  One  can 
normalize  the  probability  density  function  defined  by  these  three  values  to 
unit  variance,  for  |z(  _<  zoB  ,  by  introducing  a  fourth  parameter,  SI  .  This 
is  done  by  changing  the  variable  of  integration  from  z  to  z  as  follows: 


where 


i  =  z  Jq 


si  =  I/O  (Z0B) 


2 .  ,  if°B  2  ... 

a  (zoB)  =  J  z  pn(z)  dz 


roB  2  .  . 

=>  1  =  J  2  Pn(2)  d2  » 


JoB  -  zoB  ^ 


*  iJ5  ?. 


/  1 1® 
00  (-1J  A 

r'  a 


/ma  +  1\  |ma  +  1  1  -z2  \ 

V  2  /r  1  \  2  ;  2  ’  Q  ; 


Typical  density  functions  normalized  for  unit  power  with  zQg  =  10  are 
plotted  on  a  log-log  scale  in  figures  4  and  5.  Note  that  for  constant  Aa 
(figure  4),  the  more  impulsive  densities,  i.e.,  those  whose  tails  approach 
zero  more  slowly,  are  characterized  by  smaller  a  values.  In  fact,  one  can 
see  in  figure  5  that  the  slope  of  the  tails  is  determined  by  the  a  parameter. 
Although  the  role  of  is  not  as  readily  apparent  as  that  of  a,  it  is 
related  to  the  slope  of  the  density  function  during  the  transition  from  the 
Gaussian  portion  of  the  distribution  to  the  non-Gaussian  tails  and  to  the 
range  over  which  the  density  function  is  Gaussian-like.  The  "bendover"  point, 
z0g,  is  that  point  beyond  which  equation  (3)  no  longer  adequately  describes 
the  true  density. 
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COMPARISON  OF  NONLINEAR  DEVICES  IN  CLASS  B  NOISE  ENVIRONMENT 


The  improvement  factor  of  a  nonlinear  receiver  is  a  function  of  the  total 
noise  power,  the  shape  of  the  noise  density  function  that  yields  that  power, 
and  the  input-output  characteristics  of  the  nonlinear  device.  Equation  (4) 
indicates  that  the  noise  power  depends  on  z0u,  the  bendover  point  in 
Middleton's  Class  B  Noise  Model.  Intuitively,  the  more  zQg  exceeds  the 
suppression  range  of  a  nonlinear  device,  the  greater  the  improvement  factor 
for  that  device.  An  example  of  this  can  be  seen  in  figure  6.  Note  that  the 
increase  in  the  improvement  factor  is  entirely  due  to  the  increase  in  the 
variance  once  z0g  exceeds  a  particular  value. 

Figures  7,  8,  9,  and  10  show  the  improvement  factor  for  the  optimal  non¬ 
linear  receiver  versus  z0g  for  various  noise  statistics.  (These  results 
were  obtained  using  equations  (2a)  and  (4)  and  normalizing  as  in  equation  (5).) 
For  these  calculations  Pn(z)  for  z  >  zQg  was  assumed  to  be  zero.  As  A& 
approaches  zero  (figure  7),  the  density  function  is  Gaussian  over  a  greater 
range  of  values  and,  therefore,  the  improvement  factor  depends  strongly  on 
the  shape  of  the  tails.  If  the  tails  are  decaying  rapidly  (larger  a),  only 
modest  improvements  are  to  be  expected.  On  the  other  hand,  if  the  tails 
approach  zero  very  slowly  (small  a),  significant  processing  gains  may  be 
achieved  using  the  optimum  nonlinear  receiver  over  a  linear  one.  As  An 
increases  (figures  8,9,  and  10),  the  noise  is  Gaussian  over  a  smalxer  range 
of  values;  and  improvements  will  occur  at  smaller  values  of  zQg.  Here,  also, 
the  shape  of  the  tails  greatly  affect  the  processing  gain.  In  general,  for 
a  given  zQg,  the  improvement  factors  increase  as  increases  and  as  a 
decreases. 

One  of  the  primary  goals  of  this  report  is  to  rate  the  performance  of 
simpler  suboptimal  nonlinear  receivers.  The  improvement  factor  for  the  hard 
limiter,  the  clipper,  and  the  hole  puncher  were  evaluated  under  a  variety  of 
noise  conditions  using  Middleton's  model  and  equations  (2b) -(2d).  Typical 
results  are  shown  in  figures  11-14.  The  performance  of  the  clipper  and  the 
hole  puncher  are  functions  of  the  "clip  level,"  i.e.,  the  percentage  of  time 
the  input  is  in  the  nonlinear  portion  of  the  device’s  operating  characteris¬ 
tics.  The  hard  limiter  is  a  special  case  of  the  clipper  in  which  the  clip 
level  is  set  at  100  percent. 

Figures  11,  12,  and  13  present  the  ratio  of  the  clipper  improvement 
factor  to  optimum  improvement  factor  versus  clip  percentage  for  various  A 
and  a.  While  the  optimum  clip  level,  i.e.,  the  clip  level  that  results  in 
the  best  performance,  depends  on  the  noise  parameters  defining  the  density 
function,  the  performance  of  the  clipper  relative  to  the  optimum  nonlinearity 
is  not  a  strong  function  of  the  percentage  of  time  that  the  input  is  in  the 
nonlinear  portion  of  the  transfer  function.  In  fact,  by  setting  the  clip 
level  to  obtain  90  percent  clipping,  the  clipper  performance  will  remain 
within  about  4.5  dB  of  the  optimum  over  a  large  range  of  noise  parameters. 
When  figures  7-10  are  used  as  an  absolute  measure,  the  optimal  improvement 
factor  is  large  for  just  those  values  of  a  for  which  the  clippers'  perform¬ 
ance  decreases.  As  the  optimal  performance  drops  in  value  due  to  an  in¬ 
creasing  o,  the  clipper's  performance  increases  so  that,  relatively  speak¬ 
ing,  the  clipper  is  always  close  to  the  optimal. 
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Even  simpler  than  the  clipper  is  the  hard  limiter,  which  approaches  the 
performance  of  the  clipper  without  the  requirement  cf  adaptive  control  of 
the  clip  level.  The  more  impulsive  the  noise,  the  closer  the  hard  liiaiter's 
performance  approaches  that  of  the  optimum  clipper. 

The  relative  performance  of  the  hole  puncher  (figure  14)  is  very  sensi¬ 
tive  to  both  the  clip  level  and  the  parameters  of  the  noise  density;  but, 
when  adjusted  properly,  it  essentially  can  achieve  optimal  performance.  Thus, 
while  the  transfer  characteristic  is  easily  implemented,  logic  is  required 
to  estimate  the  noise  parameters,  and  adaptive  control  is  needed  to  maintain 
the  proper  clip  level. 

The  four  nonlinear  devices  discussed  in  the  report  form  a  hierarchy  of 
receivers  based  on  complexity  and  performance.  Although  each  application 
determines  the  exact  tradeoffs,  there  are  many  cases  in  which  the  simpler 
suboptimal  devices  are  more  cost  effective  than  the  optimal. 


ELF  RECEIVER  PERFORMANCE 


In  this  section  the  general  results  obtained  in  the  previous  sections 
are  applied  to  the  ELF  range  cf  electromagnetic  communications.  Several  digi¬ 
tal  tape  recordings  of  tvpical  high  level  ELF  noise  from  the  Saipan  area  were 
analyzed  to  yield  amplitude  probability  distributions.  These  measured  distri¬ 
butions  are  compared  with  distributions  computed  from  Middleton's  Class  B 
Noise  Model,  equation  (S),  and  the  parameters  A^j,  a,  and  zDg  are  determined. 
From  these  parameter  values  and  figures  7,  8,  and  S  the  expected  performance 
for  the  optimum  and  suboptimum  nonlinearities  can  be  determined. 


Figures  15-17  are  plots  of  the  amplitude  probability  distributions  of 
typical  high  level  ELF  noise  from  Saipan.  The  noise  was  bandpass  filtered 
to  20-150  Hz  and  each  distribution  was  measured  over  400  s.  The  ordinate 
of  these  figures  is  the  percentage  of  tine  that  the  magnitude  of  the  noise 
exceeds  the  abscissa  level;  the  reference  level  of  the  abscissa  is  the  ras 
level  of  the  noise.  The  individual  squares  on  the  plots  are  points  from 
Middleton's  model.  From  these  plots  it  is  clear  that  the  amplitude  proba¬ 
bility  distributions  and  density  functions  of  actual  ELF  noise  can  be 
described  accurately  by  the  amplitude  distributions  and  densities  resulting 
fro®  Middleton's  model.  For  these  high  levels  of  ELF  noise;  the  parameter 
values  range  over  Ajj  =  1.0  -  1.5,  a  =  1.2  -  1.4,  and  z0g  approaching  40  dB. 
An  optimally  designed  nonlinear  receiver  could  be  expected  to  perform  10-20  dB 
better  than  a  linear  receiver  for  these  parameter  ranges;  however,  the  actual 
performance  will  depend  on  the  particular  combination  of  parameter  values. 

If  a  clipper  were  utilized  and  set  to  clip  between  20-80  percent  of  the  time, 
its  performance  would  be  within  about  1  dB  of  the  optimum  nonlinearity. 


The  current  Propagation  Validation  System  (PVS)  ELF  receiver  utilizes 
a  clipper  that  adjusts  itself  so  that  the  received  signal  and  noise  is  clipped 
40  percent  of  the  time.  This  rather  simple  device,  which  does  not  necessitate 
any  complex  noise  parameter  estimation,  provides  performance  that  is 
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c oap arable  to  that  provided  with  the  optima  nonlinearity.  A  detailed 
analysis  of  the  nonlinear  noise  processing,  including  the  aspects  leading 
to  the  design  and  measured  performance  utilized  in  the  Navy’s  ELF 
receivers,  is  contained  in  references  1  and  10. 


CONCLUSIONS 


The  optimum  receiver  for  detecting  known  threshold  signals  in  additive, 
white,  non-Gaussian  noise  has  been  described  and  its  performance  is  defined 
as  a  function  of  the  noise  aaplitude  probability  density  function.  The 
optimum  receiver's  performance  has  been  calculated  and  plotted  versus  various 
ranges  of  the  parameters  defining  Middleton’s  Class  B  Noise  Model.  From 
these  plots  it  is  apparent  that  the  optimum  nonlinear  receiver  can.  yield  very 
significant  performance  improvements  relative  to  the  receiver  that  is  optimum 
in  Gaussian  noise.  Implementation  of  the  optimum  nonlinearity,  however,  can 
be  rather  complex,  requiring  real-time  estimation  of  the  noise  parameters. 

The  performances  of  several  suboptimum  nonlinearities  were  calculated 
and  their  performance  relative  to  the  optimum  was  plotted.  A  properly 
adjusted  hole  puncher  yields  performance  within  1  dB  of  the  optimum;  however, 
its’  relative  performance  is  very  sensitive  to  the  percentage  of  time  the 
input  is  suppressed.  Proper  adjustment  requires  estimating  the  noise  param¬ 
eters,  as  is  also  required  for  the  optimum  nonlinearity.  The  performances  of 
the  clipper  depend  on  the  noise  parameter  values  and  the  percentage  of  clip, 
as  does  the  hole  puncher;  however,  it  is  less  sensitive  to  these  values  than 
is  the  hole  puncher.  A  simple  clipper  adjusted  to  clip  90  percent  of  the  time 
yields  performance  within  4.5  dB  of  the  optimum  over  a  wide  range  of  noise 
parameters;  the  relative  performance  improves  as  the  noise  becomes  less 
impulsive.  The  hard  limiter,  requiring  no  adjustment,  yields  performance  within 
2  g8  of  the  optimum  clipper;  the  performance  relative  to  the  optimum  clipper 
improves  as  the  noise  becoaes  more  impulsive. 

The  instantaneous  amplitude  probability  density  function  resulting  from 
Middleton's  model  fits  measured  high  level  ELF  noise  quite  closely  with 
Aq  =  1.0  -  1.5,  a  =  1.2  -  1.4,  and  zQg  «40  dB.  For  this  range  of  parameters, 
a  clipper  adjusted  to  clip  40  percent  of  the  time,  as  does  the  PVS  ELF 
receiver,  provides  performance  within  0.5  dB  of  the  optimum. 

Recommendations  for  future  work  include 

1.  Determine  optimum  receiver's  performance  sensitivity  to  inaccuracies 
in  parameter  estimation. 

2.  Extend  results  to  determine  performance  of  these  threshold 
receivers  for  nonthreshold  signals. 

3.  Compare  predicted  performance  of  the  optimum  receiver,  as  defined 
here,  at  very  low  frequency  (VLF)  with  measured  performance  of  existing  VLF 
receivers. 
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Appendix  A 


DERIVATION  OF  THE  IMPROVEMENT  FACTORS 
FOR  SEVERAL  NONLINEARITIES 


In  this  appendix  the  improvement  factors  for  the  optimum  nonlinearity, 
the  hard  limiter,  the  clipper,  and  the  hole  puncher  are  derived.  The  improve¬ 
ment  factor  is  defined  as  the  signal  to  noise  ratio  (SNR)  at  the  correlator 
output,  with  the  appropriate  nonlinearity  in  the  signal  path,  divided  by  the 
SNR  with  the  nonlinearity  removed.  The  improvement  factor  is  an  indication 
of  the  performance  improvement  or  degradation  that  is  expected  owing  to 
inclusion  of  the  nonlinear  device  in  the  signal  path.  The  test  statistic  at 
the  correlator  output  is  the  sum  of  independent,  finite  variance  random 
variables;  and  by  the  Central  Limit  Theorem  it  is  asymptotically  Gaussian 
distributed.  Hence,  the  output  SNR  is  sufficient  for  describing  receiver 
performance  for  large  N. 

The  receiver  structure  of  interest  is  shown  in  figure  A-l.  Each  input 
sample,  x*,  is  the  sum  of  a  completely  known  signal  sample,  S^,  and  an  inde¬ 
pendent  noise  sample,  nj.  The  signal  samples  are  considered  to  have  come 
from  one  of  two  binary  equal-energy  signals;  the  noise  samples  are  identically 
distributed  zero  mean  with  variance  and  have  a  symmetric  first  order 
amplitude  probability  density  function.  The  y^  are  the  outputs  of  a  zero 
memory  nonlinearity,  f(x^).  The  test  statistic,  6,  is  the  sum  over  i  of  N 
products  of  the  form  S^y*  The  SNR  at  the  correlator  output  is  defined  as 
the  ratio  of  the  square  of  the  expected  value  of  6  and  the  variance  of  5. 
Remembering  that  the  Sj  are  completely  known,  the  SNR  at  the  correlator  output 
with  the  nonlinearity  in  the  circuit  is  easily  found  to  be  given  by 

/  N  \2 
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Figure  A-l.  Optimum  Receiver  Structure  for  Threshold  Detection 


With  the  nonlinearity  removed,  the  SNR  at  the  correlator  output  is  given  by 


SNRl  = 


N  ? 

?  » 


(A-2) 


resulting  in  an  improvement  factor,  I,  of 


(A-3) 


In  order  to  evaluate  the  improvement  factor  for  the  various  nonlinear¬ 
ities  of  interest,  the  mean  and  variance  of  the  samples  at  the  nonlinear  out¬ 
put  are  n'eeded. 

For  threshold  signals  in  independent,  identically  distributed  noise,  the 
optimum  nonlinearity  takes  the  form 


yi  =  f(V  =  ~  dF  ln  Pn(V  *  (A_4D 

i 

where  p  (x-)  is  the  amplitude  probability  density  function  of  the  noise-only 
samples.  The  expected  value  of  y.  is  given  by 

7i  ■/  f-  at:  p«lxi3 

It  is  possible  to  expand  Pj,(xi  -  S^)  in  a  Taylor  series  about  xj,  and  the 
terms  in  Sj  of  order  two  and  greater  can  be  dropped  for  threshold  signals. 
Then,  using  the  fact  that  PnCx^)  is  symmetrical,  yi  is  found  to  be 


p  (x.  - 
rn  i 


S.)  dx. 
x  1 


U-5) 


y.  =  S.L 

•'l  l 


(A-6) 


L  =/_lpn’(li)]  /  pn(IiJ  • 

Using  a  similar  approach,  we  find  the  variance  of  the  nonlineur  output  to  be 

Var(y.)  =  L(1  -  S?  L)  .  (A-7) 

The  improvement  factor  for  the  optimum  nonlinearity  is,  then,  given  by 


(A-8) 


A-2 
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with  equality  as  S.  L  approaches  zero. 

The  hard  limiter  is  defined  by  a  transfer  function  such  that  y.  =  +J 
x.  0  and  y^  =  -1  for  x.  <  0.  The  expected  value  of  y.  is  given  by 


f.  -  2  I f  p  (x.  -  S.)  dx.  - 
l  I J  rnv  l  ix 

L  ° 


(A-9) 


where  the  expression  within  the  brackets  is  simply  the  probability  of  y.  =  -*-1 
Since  we  are  interested  in  comparing  the  various  non linearities'  under  similar 
signal  conditions,  i.e.,  threshold  signals,  Sj  can  be  assumed  small.  For 
small  S^  and  a  symmetrical  noise  density  function. 


yi  =  2Si  ?n(0> 


(A- 10) 


Vax(y.)  =  1  -  4S2  p2(0)  . 


CA-11) 


Substituting  equations  (A-10)  and  (A-ll)  into  (A-3)  yields  the  improvement 
for  the  hard  limiter. 


(A-12) 


with  equality  as  4S.  Pn(0)  approaches  zero. 

The  clipper  is  described  by  a  transfer  function  such  that 


yi  =  +C 


-  C  <  X.  <  c 
—  1  — 


X.  >  C 
X 


xi  <-c* 


(A-13) 


where  c  is  a  constant  greater  than  zero.  The  expected  value  of  y.  is  found 
from  1 


y,  =  f  p_(x,  -  sj  dx,  +  f 


c  p  (x.  -  S.)  dx. 


f-c)  v  (x.  -  S.)  dx 
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Again,  letting  become  small. 


y.  =  S.p  , 
1 


where 


P  =  2  J  Pn(V  toi 


(A-15) 


In  a  sinilar  manner,  the  variance  of  the  clipper  output  can  be  found: 


Var(y.)  =  2  J  x2  pn(x.)  dx.  +  c2(l  -  p) 


(A-16) 


Substituting  the  above  expressions  into  (A-3)  yields  the  improvement  factor 
for  the  clipper: 


2  f  x2  p  (x.)  dx.  +  cz{l  -  p) 
J  i  rn  l  l  r 


(A-17) 


The  hole  puncher  is  similar  to  the  clipper  except  that  the  output  is  set 
to  zero  if  the  input  exceeds  the  threshold  c.  The  mean  and  variance  at  the 
output  of  the  hole  puncher  are  found  by  following  the  procedure  utilized  with 
the  clipper.  The  mean  is  given  by 


n  *  si  [p  -  2c  p»Cc)]  • 


(A-18) 


where  p  was  defined  in  equation  (A-15)  and  p,j(c)  is  the  noise  density  function 
evaluated  at  c.  The  variance  is  given  by 


r c  2 

Var(y.)  =  2  J  x.  Pn(x.)  dx. 


and  the  improvement  factor  by 


THP 


.  °n  [p  -  2c  Pn(c)]: 

2  f  x?  p  (x.)  dx. 
J  i  rn  x  i 


(A-19) 


(A-20) 


A-4 


It  is  interesting  to  note  that  while  the  SNR's  at  the  correlator  output 
depend  on  the  signal  structure, the  improvement  factors  for  the  four  nonlin¬ 
earities  of  interest  are  independent  o.T  the  values  of  the  signal  samples. 
This  is  a  result  of  the  small  SNR  assumption. 
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Appendix  B 


NUMERICAL  COMPUTATION’  OF  THE 
HYPERGEOMETRIC  AND  GAJMA  FUNCTIONS 


Inherent  within  the  results  of  this  report  is  the  need  to  evaluate  the 
density  function  ,  pn(z),  for  Middleton’s  Class  B  Noise  Model  and  its  deriv¬ 
ative,  p^  (z).  These  functions  are  expressed  in  equations  (B-l)  and  (B-2); 
p^  (z)  is  derived  from  pn(zj  with  the  aid  of  the  following  identities: 


£  lFi<a;b;z)  =  f-  ^(a  +  1;  b  +  1;  z) 
aF(a)  =  T(a  +  1) 
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where 


o  <  a  <  2  A>0 

a 

T(a)  =  Gamma  Function 

jF, (a;  b;  c)  =  Confluent  Hypergeometric  Function. 


Thus,  both  p  (z)  and  its  derivative  require  computating  the  Hypergeometric 
Function,  jFi(a;  b;  -z),  and  the  Gamma  Function,  T(a).  For  all  values  of  a 
and  the  index  of  summation,  a,  indicated  in  equations  (B-l)  and  (B-2),  the 
parameters  a,  b,  and  z  satisfy  a  ^  b  >  0.  and  z  >  0.  This  appendix  outlines 
the  algorithms  used  in  this  report  to  compute  the  Hypergeometric  and  Gamma 
Functions  for  these  ranges  of  arguments. 

The  Confluent  Hypergeometric  Function,  .Fjfe;  b;  -z),  is  defined  in 
equation  (B-3).  This  equation,  however,  is  suitable  for  the  numerical  com¬ 
putation  of  iFj(a;h;  -z)  only  for  small  arguments  of  a  and  z.  Thus,  equa¬ 
tions  (B-4)  through  (B-6)  are  needed  to  evaluate  the  Hypergeometric  Function 
for  a  broader  range  of  arguments,  in  particular  large  Zi 
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*  (a)  (-1)1  z1 

iV=*  “z)  =£  -  CB-3) 

1=0  (b) .  1! 

i 

jFjCajb;  -z)  =  e  -2  jF^b  -  a;  b;  z)  (B-4) 


(b  -  a)  jFjfa  -  1;  b;  -z)  +  (2a  -  b  -  z)  ^(a;  b;  -z)  -  a  ^(a  +  l:  b;  -z)  =  0 

(B-5) 

N 

iaj.  ia  +  4  -  qj- 

(2  »  0), 


F  fa*  b*  -zl  *3?  2~a  V* 

rr3' ^  J  r(b  -  a)  ** 


U)i  ta  *  *  -  »i  _ _ (B-6) 


i=o  i!  z 


where 


i-1  A 

(a).  =  a  (a  +  1)  (a  +  2)  • -  (a  +  i  -  1)  =  ir  (a  +  j);  (a)Q  =  1 

j=o 


N  =  index  of  the  smallest t era  in  equation  (B-6). 


Whereas  equations  (B-3)  through  (B-5)  are  identities,  equation  (B-6) 
is  an  asymptotic  series  approximation  to  the  Hypergeoaetric  Function  valid 
only  for  large  z  with  respect  to  a  and  b.  Note,  however,  that  this  asymptotic 
series  is  divergent  for  any  given  z  since 


lim  ,  Si+l 
n-**>  '  t 

n 


lim  Si+l  _  lim 

n-**>  t  ~  n+® 
n 


(Wntl  ta  ^  1  -  b)n<  ' 
\  (n  +  x);  2n<1  / 

i  ca)n  ta  *  1  -  bU 


lim  (a  +  n)  (a  +  1  -  b  +  n) 
n+«  (n  +  1)  z 


Hence,  the  point  of  truncation  of  the  series,  N,  is  important  when  using 
equation  (B-6)  to  approximate  the  Hypergeometric  Function. 

To  understand  the  behavior  of  the  terms  of  this  asymptotic  expansion, 
consider  the  following: 
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Since  (a)n(a  +  1  -  b)  /(n  +1)!  zn+*  is  positive,  the  sign 
determined  by  the  quadratic  equation 

[(a  +  n)  (a  +  1  -  b  +  n)  -  (n  +  1)  z] 

2  2 
=  [n  +  (2a  +  1  -  b  -  z)n  +  (a  +  a  -  ab  + 

which  has  roots  at 

Vn2  =  tz  +  b  -  2a  -  1  +  V?  +  2z(l  +  b  -  2a)  ♦  (1  -  b)2)/2  . 

Figure  (B-l)  is  a  graphical  representation  of  this  quadratic  function. 
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Figure  B-l.  Graphical  Representation  of  the  Quadratic 
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If  nj  and  03  are  both  less  than  zero,  equation  (B-6)  diverges  immediately 
and  should  not  be  used  to  approximate  the  Hypergeometric  Function  for  the 
given  values  of  a,  b,  and  z.  If,  however,  nx  and  n2  are  of  the  opposite  sign, 
the  terms  of  the  series  initially  decrease  in  value  until  n  =  n2  and,  then, 
grow  indefinitely.  When  both  nx  and  n2  are  positive,  the  terms  Initially 
increase  in  value,  then  decrease,  and  finally  grow  indefinitely  again.  In 
these  latter  two  cases,  one  can  find  the  smallest  term  of  the  series  (i.e.,  the 
stopping  point)  by  comparing  the  (n  +  l)tb  term  with  the  ntb  term  only  after 
the  inde~  exceeds  nx-  Thus,  if  n  >  nx,  tn+x  >  tn  and  tn  <  t„  =  1  are  all 
satisfied,  the  stopping  index  N  equals  n.  If  this  last  condition,  tn  <  t  ,  is 
the  only  one  not  satisfied,  the  first  term  of  the  expansion,  tD,  must  be  the 
smallest  and  N  is  set  equal  to  0. 


The  computation  of  pj,(z)  or  p^(z)  for  a  particular  z  requires  evaluating 
the  Hypergeometric  Function  via  a  series  in  which  the  "a"  argument  ((ma/2) 

+  b)]  is  a  linear  function  of  the  index  of  summation  and  is,  therefore, 
increasing  as  more  terms  are  included  in  the  series.  For  large  z,  one  hopes 
that  the  asymptotic  expression  (equation  (B-6))  can  be  used,  and  that  the 
A®  F(*)/m!  factor  in  equations  (B-l)  and  (B-2)  decreases  rapidly  enough  so 
that  only  a  few  terms  in  the  suasation  need  be  evaluated  to  satisfactorily 
approximate  Pn(z)  or  pcf(z).  By  choosing  to  use  equation  (B-6)  only  when  z  >  10, 
one  can  guarantee  that  the  asymptotic  expression  will  be  valid  to  approximate 
lFx(a;b;  -z),  for  m  =  0,1,2,  a  =  (ma/2)  +  b,  b  =  1/2  or  3/2,  and  0  <  a  <  2. 

If,  however,  the  A*  f (♦)/*!  factor  is  not  decreasing  rapidly  enough  to  limit 
the  nuabei  of  terms  required  in  equation  (B-l)  or  (B-2) ,  so  that  a  =  [  (ma/2)  +  b] 
becomes  too  large  to  use  equation  (B-6)  to  approximate  the  Hypergeometric 
Function  for  the  given  z,  then  equation  (6-6)  may  be  used  on  a  •  1  or  a  -  2  — , 
and  then  equation  (B-5)  can  be  utilized  to  iterate  on  "a”  back  to  its  original 
value. 


B-3 


TR  6339 


This  iteration  technique  was  also  used  to  evaluate  the  Gamma  Fraction, 
r(z),  which  appears  with  positive  arguments  in  equations  (B-l)  and  (B-2)  and 
negative  arguments  in  equation  (B-6).  A  polynomial  approximation  to  T(z) 
(equation  (B-7))  was  used  on  the  fractional  part  of  the  argument,  thereby 
requiring  one  to  use  the  recursion  relationship  in  equation  (B-8)  to  step  the 
argument  up  or  down  to  its  desired  value: 

8 

T(z  +  1)  ~  1  +T.  c.  z1  0  <  z  <  1  (B-7) 

i=l  1 


zr(z)  =  r(z  +  l) 


(B-8) 


r(z)  =  e~Z  z(z-3^  ^2?  e(1/12z)  30  <  z  £  55  ,  (B-9) 

where 


r> 

II 

-0.577191652 

c5  =  -0.75670  4078 

C2  = 

0.988205891 

c6  =  0.48219  9:  ; 

C3  = 

-0.897056937 

c,  =  -0.19352  7818 

n 

■o. 

ii 

0.918206857 

cg  =*  0.03586  8343 

Figures  B-2  and  B-3  are  flow  diagrams  of  the  algorithms  used  in  this 
report  to  evaluate  the  Hypergeometric  and  the  Gamma  Fractions,  respectively. 
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Figure  B-3.  Flow  Diagram  for  Evaluating  the  Gamma  Function 
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